function Shell_mixed_elecavg_grID_wtp_tomlab_varagin_j_bs(subsample,inc,set_seed,matlab_seed,get_nb_draws)

%matlab_seed=1353;
rng(matlab_seed)

	tmp_file1=strcat('Z:\FAC\HEC\DEEP\shoude\default\D2c\eegap\EEgap_data_code_heter_SM\Data\Matlab_estimation\choiceset_identifier_trimester_week_zipcode_2008_2012_v11022017_struct_',num2str(subsample),'_',num2str(inc),'_seed_',num2str(set_seed),'.csv');
    id=load(tmp_file1); 
	Nmax=size(id,1);
	bs_max=10*ceil(Nmax/20000);

index=[13,14,16];
bs_max=length(index);

parfor i=1:bs_max
    j=index(i);
	%index_bs=[(1+(j-1)*(1000)):min(j*(1000),Nmax)];	
	index_bs=[(1+(j-1)*(2000)):min(j*(2000),Nmax)];
	Master_hom_elecavg_grID_wtp_tomlab_varagin_j_bs(subsample, inc, set_seed, j, index_bs);
	Master_mixed_elecavg_grID_wtp_tomlab_varagin_j_bs(subsample,inc,set_seed,matlab_seed,get_nb_draws, j, index_bs);
end

end